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Abstract 

The forward photoproduction amplitude for Tp is calculated in a pQCD A;_L-factorization 

approach with an unintegrated gluon distribution constrained by inclusive deep-inelastic structure 
functions. The total cross section for diffractive Ts is compared with a recent HERA data. We 
also discuss the 2S/\S ratio in diffractive T-production. The amplitude is used to predict the 
cross section for exclusive T production in hadronic reactions. Differential distributions for the 
exclusive pp pTp process are calculated for Tevatron energies. We also show predictions for 
LHC. Absorption effects are included. 
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I. INTRODUCTION 



The inclusive production of quarkonia was studied intensively in the past both in ele- 
mentary hadronic and nuclear reactions at SPS, RHIC and Tevatron energies. For a review 
see e.g. In contrast, the exclusive production of heavy QQ vector quarkonium states 
(e.g. hih2 — > h{Th2) in hadronic interactions was never measured, but attracted recently 
much attention from the theoretical side [2I, Is], 0, 0, [sl, 0, [sj . Due to the negative charge- 
parity of the vector meson, the purely hadronic Pomeron-Pomeron fusion mechanism of 
exclusive meson production is not available, and instead the production will proceed via 
photon-Pomeron fusion. A possible purely hadronic mechanism would involve the elusive 
Odderon exchange [2I, [6|. Currently there is no compelling evidence for the Odderon, and 
here we restrict ourselves to the photon-exchange mechanism, which exists without doubt, 
and must furthermore dominate any hadronic exchange at very small momentum transfers. 
In our approach to the exclusive hadronic reaction, we follow closely the procedure outlined 
in our previous work on J/ip production 0]. There is one crucial difference, though. While 
in the case of diffractive J / ip photoproduction there exist a large body of fairly detailed data, 
including e.g. transverse momentum distributions, the photoproduction data for exclusive 
T's are rather sparse jo], [l^, 11|. Hence, different from 0] we cannot avoid modelling the 



relevant 7p — > Tp subprocess. Fortunately, due to the large mass of the T's constituents, 
the cross section gets its main contribution from small-size 66-dipoles, and the production 
mechanism can be described in a pQCD framework (for a recent review and references, see 
(1^). The two main ingredients are the unintegrated gluon distribution of the proton, and 
the light-cone wave function of the vector meson. The unintegrated gluon distribution is 
sufficiently well constrained by the precise small-x data for the inclusive proton structure 
function, and we shall content ourselves here with a particular parametrisation which pro- 



vides a good description of inclusive deep inelastic scattering data [13j. As the relevant 
energy range of the 7p Tp subprocess at Tevatron overlaps well with the HERA energy 
range, any glue which fulfills the stringent constraints of the precise HERA F2-data must 
do a similar job. Alternative unintegrated gluon distributions are discussed for example in 



The current experimental analyses at the Tevatron [15[ call for an evaluation of differential 
distributions including the effects of absorptive corrections. 

The HERA data cover the 7p center of mass (cm-) energy range W ~ 100-^200 GeV. This 
energy range is in fact very much relevant to the exclusive production at Tevatron energies 
for not too large rapidities of the meson, say \y\ < 3. This will be different at the LHC, 
where a broad range of subsystem energies W-yp, up to several TeV, is spanned for T emitted 
in the forward directions. This will require a long-range extrapolation to a completely new 
unexplored region. In this paper, however, we will concentrate on predictitions for Tevatron 
energies. Here our input amplitude is constrained by the HERA data, to which description 
we now turn. 



II. PHOTOPRODUCTION -/p Tp AT HERA 

We thus turn to the analysis of the photoproduction recation studied at HERA. The 
photoproduction amplitude will then be the major building block for our prediction of 
exclusive T production at the Tevatron. 
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A. Amplitude for 'jp — > Tp 




FIG. 1: A sketch of the exclusive "yp Tp amplitude. 

The amplitude for the reaction under consideration is shown schematically in FiglH As 
it is explained in Ref. the imaginary part of the amplitude for the ■y*p Tp process 
can be written as 



47r2 
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dz(Pk 
z(l-z) 



[z,k,Ki,K2,Q^) , (2.1) 



where the transverse momenta of gluons coupled to the QQ pair can be written as 



1^2 



K 



(2.2) 



The quantity jF(xi, X2, a^i, ^2) is the off diagonal unintegrated gluon distribution. Explicit 
expressions for I\^^\y can be found in For heavy vector mesons, helicity-flip transi- 

tions may be neglected, and we concentrate on the s-channel helicity conserving amplitude, 
= Xv- In the forward scattering limit, i.e. for A = 0, azimuthal integrations can be 
performed analytically, and we obtain the following representation for the imaginary part of 
the amplitude for forward photoproduction 7p Tp : 



^mMiW,A^ = 0,g2 = 0) = |^2CTV47raem 
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and 



+ ml ^(A;2-m^-K2)2 + 4m^F 
A;2 + ml f ^ k"^ — ml — 



Wi(e,K^) = 1- '" (1+ . =)• (2.7) 

To obtain these results, the perturbative 7 ^ 66 hght cone wave function was used; the 
vertex for the bb ^ T transition is given below, and is obtained by projecting onto the 
pure s-wave 66-state. Here cx = Cf, = — 1/3, and the mass of the bottom quark is taken as 
rrib = 4.75 GeV. The relative transverse momentum squared of (anti-)quarks in the bound 
state is denoted by k"^, their longitudinal momentum fractions are z,l — z, and we introduced 

= . (2.8) 

The dominant contribution to the amplitude comes from the piece cx AqVFo ~ ^IWo, and 
our exact projection onto s-wave states differs in fact only marginally from the naive 7^- 
vertex for the T ^ bb transition. The 'radial' light-cone wave function of the vector meson, 
ipvizjk"^) will be discussed further below. The unintegrated gluon distribution J-'{x,k,'^) is 
normalized such that for a large scale it will be related to the integrated gluon distribution 
g{x, Q"^) through 

xgix,Q') = ^J^(x,k'). (2.9) 

The running coupling as enters at the largest relevant virtuality g2 = max{K2, A;2 + ml}. 
Due to the finite mass of the final state vector meson, the longitudinal momentum transfer 
is nonvanishing, and, as indicated above, a more precise treatment would require the use of 
skewed/off-diagonal gluon distributions. At the high energies relevant here it is admissible 
to account for skewedness by an appropriate rescaling of the diagonal gluon distribution [13] . 



With the specific gluon distribution used by us, the prescription of [17j can be emulated by 
taking the ordinary gluon distribution at l2| 



M2 M2 

Xe// = ~ 0.41 ■ (2.10) 

The full amplitude, at finite momentum transfer, well within the diffraction cone, is finally 
written as 

M{W, A^) = {i + p) ^mM{W, = 0) exp{-B{W)A'^) . (2.11) 

Here A2 is the (transverse) momentum transfer squared, B{W) is the energy-dependent 
slope parameter: 

BiW) = Bo + 2a',jf\og[^), (2.12) 

with a'^jj- = 0.164 GeV^2 [X5j^ = 95 GeV. For the value of Bq see the discussion of the 
numerical results below. For the small size bb dipoles relevant to our problem, a fast rise of 
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the cross section can be anticipated, and it is important to include the real part, which we 
do by means of the analyticity relation 
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Finally, our amplitude is normalized such, that the differential cross section for — >■ Vp 
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B. bb wave function of the T meson 

We treat the T, T' mesons as bb s-wave states, the relevant formalism of light-cone 



wavefunctions is reviewed in 12 . The vertex for the T — > 66 transition is taken as 



e, u{p,)T^v{p-,) = [M' - M^] M^, k') u{p,) (7^ - ^\^^) ^iPl) , (2-16) 

where is the polarization vector of the vector meson V = T, T'. and p'^j^ are the on- 
shell four-momenta of the b,b quarks, p\i = iTLf. The so-defined radial wave-function 
ipvizjk'^) can be regarded as a function not of z and k"^ independently, but rather of the 
three-momentum p of, say, the quark in the rest frame of the bb system of invariant mass 
M, p= {k, {2z - l)M/2). Then, 

, , , dzd'^k 4 d^p „ — 4m? 

^y{z,k') ^ ^y{p^) ^^^p^ = 1_ 2.17 

z[\ — z) M 4 

We assume that the Fock-space components of the T, T'-states are exhausted by the bb com- 
ponents and impose on the light-cone wave function (LCWF) the orthonormality conditions 
(z,j = T,r): 

= NcJ -0^^AMUp')^,{p') . (2.18) 



Important constraints on the LCWF are imposed by the decay width V — > e+e 

4:iTa^^c% o . 16 

where ([l3,[l8j) 



r(V ^ e+e-) = ^";;73 ^ ■ 9l ■ Knlo , K^lo = 1 - ir'^siO , (2-19) 



We f d^p , , , , 2^ 



9v = ^ I J^siM + m,) Mr) ■ (2.20) 
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For the - fully nonperturbative - LCWF we shall try two different scenarios, following again 
the suggestions in (12I . 18| . Firstly, the Gaussian, harmonic-oscillator-like wave functions: 



i^^sip') = Ci exp , ^p,s{p') = C,{^o - p'al) exp , (2.21) 

and secondly, the Coulomb-like wave functions, with a slowly decaying power-like tail: 



i^isiP = -j= I, , 2 2^2 ' ^2SU0 = -j= 7^ — ■ 2.22 
a/M (1 + alp^Y VM (1 + a^2P^-^ 



The parameters af are obtained from fitting the decay widths into e+e~, whereas ^o, and 
therefore the position of the node of the 25* wave function, is obtained from the orthogonality 
of the 25* and IS* states. We used the following values for masses and widthes: M{T{1S)) = 
9.46 GeVJW'(T(25)) = 10.023 GeV, r(T(15) ^ e+e-) = 1.34 keV, r(T(25) ^ e+e") = 
0.61 keV H. 



C. Numerical results and comparison with HERA data 




FIG. 2: atotilV — > '^{^S)p) as a function of the 7p cm-energy versus HERA-data. Left: dependence 
on the treatment of the bb ^ T transition; solid curves: Gaussian (G) wave function, dashed curves: 
Coulomb-like (C) wave function. Thick lines were obtained including the NLO-correction for the 
T decay width, while for the thin lines K^lo = 1- Right: dependence on the slope parameter 
(given in GeV~^), for the Gaussian wave function. The experimental data are taken from [s, 10, 11| 



In Figl2]we show the total cross section for the exclusive jp — > Tp process as a function of 
the 7p cm-energy. In the left panel we show results for two different wave functions discussed 
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FIG. 3: The 25'/lS'-ratio atotilV — > {'^S)p) / utotilV '^{^S)p) as a function of the 'jp cm-energy. 



in the text: Gaussian (solid hues) and Coulomb-like (dashed lines). Free parameters of the 
wave function have been adjusted to reproduce the leptonic decay width in two ways: (a) 
using leading order formula (thin lines) and (b) inlcuding QCD corrections (thick lines). 
Including the KNLO^^^ctor in the width enhances the momentum-space integral over the 
wave function (the WF at the spatial origin), and hence enhances the prediction for the 
photoproduction cross section. Notice that strictly speaking inclusion of the a^-correction 
is not really warranted given that we do not have the corresponding radiative corrections to 
the production amplitude. Fortunately, due to the large scale ml, the ambiguity in the two 
ways of adjusting the wave function parameters leads to only a marginal difference in the 
total cross section over most of the relevant energy range. To be fair, it should be mentioned, 
that the situation with the next-to-leading order corrections to diffractive vector mesons 
is not a very comfortable one, see for example the instabilities reported in [20]. But then, 
the systematic extension of /c_L-factorisation is yet lacking, so that at present we must be 
content with estimates of the theoretical uncertainties obtained by changing the principal 
parameters in the calculation. 

As can bee seen from the figure, different functional forms of the LCWF can lead to 
a quite substantial differences in the predicted cross section. Finally, the absence of ex- 
perimental data for t-distributions leaves the slope parameter Bo only badly constrained. 
The full, energy dependent slope can be decomposed into three contributions: one from the 
transition 7 — V^, a second one from the dynamics of the gluon ladder exchanged - which 
induces the main part of its energy dependence, and a third one from the elastic p ^ p 
vertex. In comparison to J/V'-production, we may expect, that the slope in our case re- 
ceives a smaller contribution from the 7 — V transition, due to the smaller transverse sizes 
involved |21]. It may therefore be expected that should be somewhat smaller than in 



J/\E' photoproduction, where it is ~ 4.6 GeV ^ [16 
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We show the sensitivity to the slope parameter Bo in the right panel of Figl2l 

We observe, that in general our predictions are systematically somewhat below the ex- 
perimental data. In principle, the agreement could be improved by choosing an abnormally 
small value for Bq, we shall however refrain from such an option. In our view the description 
of data, given the large error bars, is quite acceptable. The energy dependence of our result 
corresponds to an effective Ajp ~ 0.39. For our predictions for Tevatron we shall use the 
Gaussian LCWF option, with the NLO correction to the width included. 

In Figl3]we show the ratio of the cross section for the first radial excitation T(25') to the 
cross section for the ground state T(IS'). The principal reason behind the suppression of the 
25* state is the well-known node effect (see |22| and references therein) - a cancellation of 
strength in the 25* case due to the change of sign of the radial wave function. It is perhaps 
not surprising, that the numerical value of the 25'/lS'-ratio is strongly sensitive to the shape 
of the radial light-cone wave function. 

Here we assumed an equality of the slopes for T(15') and T(2S') production. This appears 
to be justified, given the large spread of predictions from different wave functions. We finally 
note, that the ratio depends very little on the choice of the Kjsilo factor (compare left and 
right panel). 

III. EXCLUSIVE PHOTOPRODUCTION IN pp COLLISIONS 
A. The absorbed 2^3 amplitude 

The necessary formalism for the calculation of amplitudes and cross-sections was outlined 
in sufficient detail in Ref. 0]. Here we give only a brief summary. The basic mechanisms 
are shown in FigHl The major difference from HERA, where the photon was emitted 




FIG. 4: A sketch of the two mechanisms considered in the present paper: photon-pomeron (left) 
and pomeron-photon (right), including absorptive corrections. 

by a lepton which does not participate in the strong interactions, now, both initial state 
hadrons can be the source of the photon. Therefore, it is now necessary to take account 
of the interference between two amplitudes. The photon exchange parts of the amplitude, 
involve only very small, predominantly transverse momentum transfers. In fact, here we 
concentrate on the kinematic domain, where the outgoing protons lose only tiny fractions 
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Zi,Z2 <^ 1 of their longitudinal momenta, in practice z < 0.1 means y < 3. In terms of 
the transverse momenta of outgoing hadrons, Pi 2, the relevant four-momentum transfers 
are = —{p^ + zfm'^)/{l — Zi) ,i = 1,2, and si ~ (1 — Z2)s and S2 ~ (1 — zi)s are 
the familiar Mandelstam variables for the appropriate subsystems. Photon virtualities 
are small (what counts here is that <C M^), so that the contribution from longitudinal 
photons can be safely neglected. Also, as mentioned above, we assume the s-channel-helicity 
conservation in the 7* — > T transition. In summary we present the 2 — > 3 Born-amplitude 
(without absorptive corrections) in the form of a two-dimensional vector (corresponding to 
the two transverse (linear) polarizations of the final state vector meson): 

^^°nPi,P2) = ei-^J'x[xAPi,h)M^*H,^vh,{s2,t2,Qi) + (1^2) (3.1) 

Zi ti 

Inclusion of absorptive corrections (the 'elastic rescattering') leads in momentum space to 
the full, absorbed amplitude 

M{p,,P2) = J 0^^S,i{k)M^'\p,-k,p, + k)=M'^'\p,,p,)-6M{p,,p,). 

(3.2) 

With 

S^iik) = {2n)'5^'\k) - \T{k) , T{k) = aZ{s) exp ( - ^B^^k') , (3.3) 
where crff((s) = 76 mb, B^i = 17 GeV~^ jl^ , the absorptive correction 6M reads 

6Mip„p2)= J ^|^T(fc)MW(pi-fc,P2 + fc). (3.4) 
The differential cross section is given in terms of M as 

da = -^{Ml"^ dydtidt2d(f) , (3.5) 

where y is the rapidity of the vector meson, and is the angle between p^ and Pg- 



B. Results for Tevatron 

We now come to the results of differential cross sections for T production. In FigO we 
show the distribution in rapidity of T(IS') (left panel) and T{2S) (right panel). The ratio 
between 25* and IS* follows closely the photoproduction ratio discussed in Sec. Ill CI The 
parameters chosen for this calculation correspond to the Gaussian wave function, with K^lq 
included in the adjustment to the decay width. Also the unintegrated gluon distribution 
is the same as the one used in section III CI The results obtained with bare amplitudes are 
shown by the thin (red) lines, and the results with absorption effects included are shown by 
thick (black) lines. Here the absorption effects are truly a correction and cause only about 
20-30% decrease of the cross section. This is in sharp contrast to the situation for the fusion 
of two QCD ladders (relevant for the production of scalar charmonia or Higgs boson). The 
rapidity distribution is only slightly distorted by absorptive corrections. Notice that larger 
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FIG. 5: Differential cross section da/dy for T(IS') (left panel) and T{2S) (right panel) for the 
Tevatron energy W = 1960 GeV. The thin solid line is for the calculation with bare amplitude, 
the thick line for the calculation with absorption effects included. 

rapidities mean also larger photon virtualities and therefore somewhat smaller transverse 
distances in the pp collision are relevant. 

Finally, in the following figures we show distributions of T's in transverse momentum. We 
show results for different values of rapidity: y = (solid), y = 2 (dashed) and ?/ = 4 (dotted). 
In Figl6] we show the distributions for T(IS') and in Fig|7| for T{2S). Both, results with 
bare amplitudes (left panels), and with absorption (right panels) are shown. The inspection 
of the figures shows that absorption effects are larger for large values of the T transverse 
momenta - they can lower the cross section by almost an order of magnitude at the largest 
transverse momenta. There is again a different effect of absorption for different rapidities. 

Notice, that our predictions, which use the \ow-z approximation of the photon flux are 
most accurate at ?/ < 3. This is quite appropriate for Tevatron, where it seems that a 
measurement is possible only at rather low rapidities. We do not show here observables 
related to outgoing proton or/and antiproton as they cannot be studied experimentally at 
the Tevatron. There will be, however, such a possibility at the LHC 

There are important issues regarding the extrapolation to LHC energies. Firstly the 
energy of the 7p — > Tp process can vastly exceed the HERA range, and secondly the 
much increased rapidity range may increase the importance of high-mass diffraction for the 
absorptive corrections. Still, to give the reader a rough idea of the expected cross section, we 
show in Fig. [S] selected spectra at the LHC energy of ly = 14 TeV. Here, in the absorptive 
corrections, we used a Pomeron intercept of Ajp = 0.08. It is interesting to point out that 
the rise towards the maximum in the rapidity dstribution reflects the energy dependence 
of the 7p — s> Tp subprocess. Absorptive corrections in that subprocess, which we neglected 
so far can possibly alter the shape of the rapidity distribution. Since there are many other 
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FIG. 6: Invariant cross section da/dydp^ for as a function of for T(15) for W = 1960 GeV. 
The solid line: y = 0, dashed line: y = 2, dotted line: y = 4. Left panel: without absorptive 
corrections; Right panel: with absorptive corrections. 

interesting aspects at larger energies we leave a more detailed analysis for LHC for a separate 
publication. 

A brief comment on previous works is in order. In 0, 0, 0] absorptive corrections were 
not included. The equivalent photon approximation is used in 0, S], which allows only to 
obtain rapidity spectra. The form of the transverse momentum distribution suggested in 
0] is not borne out by our calculation. Cross sections da/dy obtained in 0, 0, [sl lie in the 
same ballpark as the results presented here. However the shape of the rapidity distribution 
in Isl is different from ours. 



IV. CONCLUSIONS 

We have calculated the forward amplitude for 7j9 Tp reaction within the formalism 
of fc^-factorization. In this approach the energy dependence of the process is encoded in 
the x-dependence of unintegrated gluon distributions. The latter object is constrained by 
data on inclusive deep inelastic scattering. The t-dependence for the ■yp — ^ Tp process 
involves a free parameter and is in effect parametrized. Regarding the 7 ^ T transition, 
we used different Ansatze for the bb wave functions. The results for T(IS') production 
depend only slightly on the model of the wave function, while the 2S/1S ratio shows a 
substantial sensitivity. We compared our results for the total cross section with a recent 
data from HERA. Our results are systematically somewhat lower than data, although the 
overall discrepancy is not worrysome, given the large uncertainties due to the rather poor 
experimental resolution in the meson mass. The amplitudes for the jp Tp process are 
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FIG. 7: Invariant cross section da/dydpf for as a function of pf for T(25) for W = 1960 GeV. 
The solid line: y = 0, dashed line: y = 2, dotted line: y = 4. Left panel: without absorptive 
corrections; Right panel: with absorptive corrections. 

used next to calculate the amplitude for the pp — > ppT reaction assuming the photon- 
Pomeron (Pomeron-photon) underlying dynamics. In the present approach the Pomeron is 
then described within QCD in terms of unintegrated gluon distributions. We have calculated 
several differential distributions including soft absorption effects not included so far in the 
literature. Our predictions are relevant for current experiments at the Tevatron, predictions 
were made - with qualifications - for possible future experiments at the LHC. 
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